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Abstract 

In this paper, we present statistical simulators for cascaded Rayleigh fading channels with and 
£SJ . without line-of-sight (LOS). These simulators contain two individual summations and are therefore easy 

to implement with lower complexity. Detailed statistical properties, including auto- and cross -correlations 
of the in-phase, quadrature components of the channels, envelopes, and squared envelopes, are derived. 
The time-average statistical properties and the corresponding variance are also investigated to justify 
that the proposed simulators achieve good convergence performance. Extensive Monte Carlo simulations 
are performed for various statistical properties to validate the proposed simulators. Results show that 
the simulators provide fast convergence to all desired statistical properties, including the probability 
| density function (PDF), various auto- and cross-correlations, level crossing rate (LCR), and average 

c/3 . fading duration (AFD). While various tests and measurements in dense scattering urban and forest 

environments indicate that mobile-to-mobile channels may experience cascaded Rayleigh fading, the 
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proposed statistical models can be applied to simulate the underlying channels. 



Index Terms 



Cascaded Rayleigh fading, line-of-sight, mobile-to-mobile, statistical channel model, envelope, auto- 
correlation, cross-correlation, level crossing rate, average fading duration. 



^ " I. Introduction 



Mobile-to-mobile communications have found more and more applications in systems where the 
terminals and nodes are no longer stationary HI. 121. such as mobile ad hoc wireless networks O, wireless 
local area networks @, intelligence transportation systems (21-0, and vehicular-to-vehicular systems [8 ]. 
Consequently, successful development of statistical models or simulators for mobile-to-mobile (M-M) 
channels is inevitable in order to understand the statistical properties of the underlying channels and 
effective designs for the systems. A pioneering mathematical reference model for M-M channels is found 
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Fairmount, Wichita, KS 67260. Emails : {yxibdah, yanwu. ding} @ wichita.edu. Part of the material in this paper was presented 
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in the earlier work of Akki and Haber Q and Akki ifTOll . This analysis is then extended to account for 
scattering in three dimensions by Vatalaro and Forcella flTj. A modified statistical model for suburban 
outdoor-to-indoor M-M communications channels has been discussed in fT2l . A double-summation-of- 
sinusoids statistical model is applied to simulate M-M channels, assuming omni-directional antennas and 
isotropic scattering around the transmitter (Tx) and the receiver (Rx) |[T3l . Due to the uniform scattering 
around the Tx and Rx, the statistical model in lfl"3l is also referred to as a "double -ring" statistical 
model (the single-ring statistical model is proposed for fixed-to-mobile channels lTT4l ). The double-ring 
model is further modified to generate multiple uncorrected complex faded envelopes using orthogonal in- 
phase and quadrature components [15]. The modification offers a faster convergence and lower variance 
at higher complexity. The double -ring concept is also applied to modeling M-M channels where a line-of- 
sight component exists between the Tx and Rx |[T6l (the single-ring-based fixed-to-mobile Rican channel 
models can be found in ifTTl . HH). In |[T9l , stochastic properties for M-M narrow-band channels are 
derived from a geometrical two-ring scattering model, assuming that both the transmitter and the receiver 
are surrounded by an infinite number of local scatters. The model is further extended to a three-ring 
scattering model 11201 . 

While the double-ring model characterizes a category of mobile-to-mobile channels that follow a 
Rayleigh fading distribution, various measurements in dense scattering urban and forest environments 
where signals diffract from street corners, building edges, and moving vehicles, suggest that if large 
separations between the Tx and Rx are much larger than both radii rings flT9l at the TX and Rx, then the 
mobile-mobile channels follow fading distributions that are more severe than a single Rayleigh fading. 
In lfT2l . the outdoor-to-indoor M-M radio channel measurements show that the received power follows a 
double-Rayleigh distribution rather than a single Rayleigh distribution. Similar observations are seen in 
the measurements [21], whereas the Tx is moving in within a circle of a 1 m and the RX are randomly 
chosen within a 10m by 10m square with no line of sight. The measurements in |[22l - |[24l indicate that 
severe (worse than Rayleigh) fading in vehicle-vehicle channels is observed in small cities, urban, and 
suburban areas, motorways, and highways. 

In 11251 . |[26l . the concept of a mixture of single, double, and triple Rayleigh distributions is introduced 
to represent a general M-M channel. The general multiple Rayleigh channel transfer function was initially 
suggested in |[T2l for studying dense scattering environments. Effective statistical simulators for cascaded 
or multiple Rayleigh fading M-M channels seem sparse in the literature. In this paper, we develop 
two statistical simulators for a cascaded Rayleigh fading channel with no line-of-sight (NLOS) and 
extend the models to channels with LOS. The simulators contain two individual summations; therefore, 
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they are easy to implement with lower complexity. With appropriately chosen parameters, the proposed 
models can achieve the desirable statistical properties in a small number of trials, and provide satisfactory 
convergence performance. Detailed statistical properties, including auto- and cross-correlations of the in- 
phase, quadrature components of the channels, envelopes, and squared envelopes are derived for the 
simulators. The time-average statistical properties and corresponding variance are also investigated to 
justify that the proposed simulators can achieve good convergence performance. Extensive Monte Carlo 
simulation results are provided for various statistical properties to validate the proposed simulators. 

The remainder of this paper is organized as follows: Section II provides a brief review of the general 
multiple Rayleigh channel transfer function and presents the proposed channel models without LOS and 
with LOS. Section III presents detailed statistical properties of the proposed models. Simulations results 
are provided in Section IV, and Section V concludes the paper with a summary of observations. 

II. Statistical Channel Simulators for Cascaded Rayleigh Fading Channels 
A. Transfer function of cascaded Rayleigh fading channels 

In highly dense scattering environments, the transmission path between the Tx and Rx are composed of 
many scattering paths. Fig. ?? illustrates a multiple scattering environment between two mobile terminals 
with no LOS exists between the Tx and Rx (first graph) and a LOS (second graph). We assume that the 
distance between the transmitter and receiver is much larger than both radii of the scattering rings at the 
Tx and Rx, the channel in the first graph follows cascaded Rayleigh distribution with NLOS, and that in 
the second graph cascaded Rayleigh distribution with LOS. 

A cascaded Rayleigh transfer function is suggested in |[T2l . assuming that the two groups of scatter 
at the transmitter and the receiver are stationary and located a relatively large distance apart with a 
separation distance of D 3> Rt + R r , where R t and R r are the radii of the scatters around the Tx and 
Rx, respectively. The transfer function is presented as 



where subscripts T and R are associated with Tx and Rx, respectively, A n ,$ Tn ,n = l,--- ,N, are 
the identically distributed random (i.i.d) amplitude and phase of N scatter components around the Tx, 
B m , & Rm , m = 1, • • • ,M are the i.i.d amplitude and phase of M around the Rx, G T and G R are antennas 
gains, k Tn and k Rm are unit vectors corresponding, respectively, to the direction of departure and arrival 
of the n-th scatter at Tx and m-th scatter at the Rx, and w Tn and uj Rm are Doppler spreads of the Tx and 
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Rx. Efficient statistical simulators for £[]) were left unexplored. In this section, we develop four channel 
simulators (Simulators A, B, C, and D) for the cascaded Rayleigh fading channels without and with LOS. 

B. Statistical channel simulators without LOS 

In this subsection, two simulators, namely Simulators A and B, are presented for cascaded Rayleigh 
fading without LOS in CD- 

1) Simulator A: In the statistical channel model, the complex scattering components in cascaded 
Rayleigh fading channels are given by 




where Q and P designate the number of scatters around the Tx and Rx, respectively, f\ and /2 are the 
Doppler's shift frequencies, 9 n and <I> m £ [— it, 7r), n = 1, • • • , Q, m = 1, • • ■ , P, are the phase shifts for 
each scatter from the Tx and to the Rx and they are i.i.d for all n and m. 7„ represents the angle of 
departure for the n-th scatter at the Tx, and Cm is the angle of arrival for the m-th scatter at the Rx, 
given by, respectively, 

2mr — 7r + ip _ 2rmr — it + ip 
ln = 4Q ' Cm = 2P ' (3) 

where ip, ip are independent and uniformly distributed in [— it, tt). While those angles can be chosen as 

either dependent on or independent of n and m |[T3l , |[T4ll . |[T6l , |fT8l , we choose them to be independent 

in order to reduce the complexity of the simulator. 

2) Simulator B: The angles of departures and arrivals for Simulator A, as indicated in ©, are in 
[0,7r/2) and [0, it), respectively. In fact, smaller angles of departures and arrivals in the scatters are 
expected for the channels under discussion, because the distance between Tx and Rx is relatively large with 
D 3> Rt + R r - However, our simulations suggest that if the range of angles of arrivals is chosen smaller, 
the convergence performance of Simulator A degrades, especially for the cross-correlation between the 
in-phase and quadrature components. We seek an alternative simulator with an improved performance and 
the angles of departures and arrivals are both in [0,7r/2) to achieve a better representation of scattering 
patterns in the channels between a TX and Rx with larger separations. 

We assume that Simulator B has the following form: g B (t) = (gic{t)+39is(t)) (d2c(t)+j92s(t)) , where 
gic(t) and gi s {t),i = 1,2, are the sinusoidal functions which characterize, respectively, the scattering at 
the Tx and Rx. Here, we present the steps to obtain the functions at Tx, g\ c (t) and g\ s (t), while same 
procedure applies to obtain <?2c(i) and g2 S (t) at the Rx. 



July 17, 2012 



DRAFT 



5 



Let ctk = 2kn ~^ + ' l l J ; an d Q k ^ f~ = 1,... , if Tx , denote the angle of departure and phase for the k- 
th scatter at the Tx, where tp, Ok are statistically independent and uniformly distributed in [—it, it) 
for all k, and Kj x is the number of scatters around the Tx. Next, we evaluate a summation of series 
of exponentials: J2k=\ ex P (jC^/ii cos(afc) + While the summation admits a similar expression 
as the first summation in Simulator A, we consider a special case, Kj x = AN with N > 1 
being an integer. The summation is obtained as four terms corresponding to the quarters of Kj x , 
J2k=i ex P (j( 27T fit cos(afc) + 9 k )) = J2e= fy. where ^ = En=w+i ex P (i(2vr/it cos(a n ) + <9 n )) = 

+ 8 n+ £^)), £ = 0, • • • ,3. Introducing uniform distributed random 
phases n £ [— it, 7r), n = 1, • • • , iV, where @ n is independent of ?/>, for all n and one can rewrite 
Qi, i = 1, 2, 3 as 

TV 

fil = ex P (j'(2^/it cos(a n + ~) - 9 n )) , (4) 

n=l 

TV 

^2 = ex P (j(2vr/itcos(a n +7r) - n )), (5) 

n=l 

tv 3?r 
fi 3 = ^ exp (jfaht cos(q„ + y) + 8 n )) ■ (6) 

71=1 

The range of angles of departure in f^, £ = 0, • • • , 3 is in (0, 7r/2) as specified by a n , n = 1, • • • , N (the 
range give by k = 1, ••• , K Tx is in (0,7r/2)). We assign the sinusoidal functions at the Tx, i.e. 
gi c {t) and gi s (t) equal, respectively, to the terms associated with cos(a n ) and sin(a„) in J2e=o ^> 
along with a normalizing factor. Notice Slo + ^2 = 2J2n=i cos (2vr/it cos(a n ) + # n ), and f2i + = 
2En=iCos(27r/itsin(a n )+e n ). We obtain 5lc (t) = ^En=l cos ( 27r A* cos(a n )+e n ), andg ls (t) = 



jj- J2n=i cos (27r/ii sin(a n ) + n ) for Simulator B. It is worth to mention that the choice of g\ c {t) and 



TV Z^n= 

ffis(t) at the Tx coincide with the terms in the statistical model of fixed-to-mobile cellular channels lfl4l . 

Let Pk = 2k7T ~^ +t f t and k = 1, • • • , K^ x denote the angle of arrival and phase for 

the k-th scatter at the Rx, where tp, &k are uniformly distributed in [— 7r,7r), and they are im- 
pendent for all k, and iT Rx is the number of scatters around the Rx. Evaluating the summation 
J2k=i ex P (i(2 7r /2* cos(/3fc) + for K Rx = 4M, and following similar procedure describe above, we 
obtain g 2c (t) = V?E^i cos (27r/ 2 i cos(/3 m ) + $ m ), and g 2s (t) = \fi§Y%=i c os (2irf 2 t sin(/3 m ) + 
^m) for Simulator B. The range of angles of arrivals at the Rx is in (0, 7r/2) as specified by 
/3 m , m = 1, • • • , M. In summary, Simulator B is expressed as 

9 B {t) = (gu(t) + jgis(t)){g 2c (t) + jg 2a (t)) (7) 
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where gi c {t) and gi s (t),i = 1,2 are denned as 



/~7f n n 

gum = Vir E = v V E (8) 



AT ^ v " a v y V iV 

n=l n=l 

I fK M I Ik m 

92c(t) = E B m (t),92s(t) = E A»(t), (9) 

m=l m=l 

A„(f) = cos (27r/itcos(a„) + 6» n ),C n (t) = cos (27r/i* sin(a n ) + 6 n ), n = 1, • • • , N, B m (t) = 

cos (27r/ 2 tcos(/3 m ) + $ m ), and D m (t) = cos (27r/ 2 t sin(/3 m ) + * m ),m = 1, • • • ,M, /i and / 2 are 

the Doppler's shift frequencies, 6>„, O n G [-tt, vr), n = 1, • • • , iV are the phase shifts in the ?i-scatter at 

the Tx, <3? m , G [— 7r, 7t), m = 1, • • • , M are the phase shifts in the m-scatter at the Rx and are i.i.d 

for all n and m, a n is the angle of departure of the n-th scatter at the Tx, and /3 m is the angle of arrival 

of the m-th scatter at the Rx, and they are calculated, respectively, as 

2mr — 7r + ib 2mir — ir + cp 
° n = AN ' Pm = AM ' (10) 

where tp, ip are independent and uniformly distributed in [— 7t,7t), and they are independent to all the 

phases at the Tx and Rx. It is indicated that a statistical M-M model may result in faster convergence 

rates by choosing a smaller range of angles of arrivals and departures (13). Our simulations justify the 

observation as well. Compared with Simulator A, Simulator B has faster convergence and requires fewer 

number of trials to converge to the desired statistical properties. 

C. Statistical channel simulators with LOS 

Adding the LOS component to Simulators A and B, we present two simulators for cascaded Rayleigh 
fading channels with LOS, namely Simulators C and D, as 

Simulator C: h a (t) = ^ ' J- , (11) 

y/2(l + K) 

Simulator D: h D (t) = gj,l ,tV 7l , (12) 

where K is spectral to the scatter power ratio, (j>Q is uniformly distributed in [— tt, it), is the Doppler 
frequency caused by the relative velocity, because both Tx and Rx have mobility, ^3 is the relative 
angle between the relative movement and the LOS component, and the values of and <p^ are given, 
respectively, by lfl6l as 



h 



\J (|ui| cos(0i 2 ) - \v2\) 2 + (\v\\ sin(0i 2 )) 



(13) 
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where v\ and t> 2 are, respectively, the speeds of the Tx and Rx, v% is the relative speed calculated by 
^3 = /3^> A is the wavelength of the carrier, <j>\ is the angle between the Tx and the LOS, and 0i 2 i s the 
angle between the Tx and Rx directions. Simulator D has a slightly additional complexity compared to 
Simulator C, due to the higher complexity in gs{t) than in gA(t). However, as indicated by the simulation 
results, Simulator D provides faster convergence to the desired statistical properties even for a lower value 
of K. 

III. Statistical Properties of Proposed Models 

Because Simulators A and B, g A (t) and g B (t), have the same statistical properties, as do Simulators 
C and D, h c {t) and h D (t), we focus on Simulators B and D when representing the statistical properties 
in this section. To make notations less bulky, we drop the subscripts, and use g(t) and h(t) to present, 
respectively, the simulators without and with LOS. 

A. Second-order statistics for Simulator B 

Let g c {t) = Re(<?(i)) and g s (t) = Im(g(i)) be the real (in-phase) and imaginary (quadrature) parts of 
Simulator B. The autocorrelation, the cross-correlation of the in-phase, quadrature components, and the 
autocorrelation of complex envelopes are given below. Steps for the proof are presented in Appendix |A] 

R gc9c {r) = E[g c (t + r)g c {t)] = J (27r/ir) Jo(2tt/ 2 t) (15) 

R gsgB {r) = E[g s (t + r)g s (t)] = R ga9c {r) (16) 

R gcgs (T)=E[g c (t + T)g s (t)] =0 1 Rg s g c (r)=K[g s (t + T)g c (t)] =0 (17) 
Rgg(T) = ^E[g(t + r)g*(t)] 

= + R 9 o 9 s{r) + Rg. go (r) + Rg s g s (r)) = M^fir) J (27r/ 2 r) (18) 



R\ g? \ g? {r) = E[g 2 c (t)g 2 c (t + r) + g 2 s {t)g 2 c {t + r) + g 2 (t)s 2 s (t + r) + g 2 (t)g 2 (t + r)] 
= 4 + 4J 2 (2, h r) + U 2 (2nf 2 r) + U 2 (2n fl r)J 2 (2n f 2 r) + M^hr)M^hr) 

Jo(47t/it) + J (47r/ 1 r)J 2 (27r/ 2 r) J (47r/ 2 r) + J (47r/ 2 r)J 2 (27r/ 1 r) 4g(/ 1; r)g(/ 2 , r) 

N M N 2 M 2 
4M(l + J 2 (27r/ 2 r)) + Jo(4vr/ 2 r) 4iV(l + J^hr)) + J (47r/ir) 
WM e(/l ' T) NAP e(/2 ' T) (19) 
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where Jo is the zero-order Bessel function first kind, £(/i,t) = J2n=i (e[cos (27r/ir cos(q„))] j , 
and S,(/2,t) = J2m=i (e[ cos (2vr/2r cos(/3 m ))] J . The expectation terms in (fl9l ) can be evaluated 



numerically. It is also worth noting that, although Simulator B has a different probability density 
distribution function than the double -ring simulator in 021, their autocorrelation and cross-correlation 
properties are the same. 

The time-average correlations for Simulator B can be derived as 

1 /•' ,.. .. ... 1_ 

n,m=l 



^ (t)= t™t/ Sc(t)g c (t + r)dt = ^^ J2 {An(r)Bm(T) + C n (T)D m (r)) (20) 

n i-n — 1 



1 



T 



^W T ) = J im 7f / 9c(t)g s (t + r)dt = (21) 







Rg.g.(T) = Jim ^ ^ 9s{t)g c (t + r)dt = (22) 







1 - T N,M 

R 9sgs (r) = lim^- J g s (t)g s (t + r)dt= mw ]T (A»(r)D TO (r) + C n (r)B m (r)) (23) 



n,m=l 



10 

where A n {r) = cos (27r/ir cos(a n )), -B m (r) = cos (27r/ 2 r cos(/3 m )), C n (r) = cos (27r/iTsin(a n )J, and 

D m {r) = cos (27r/2r sin(/3 m )). The variance of the time-average correlations for Simulator B can be 
derived as 

Var[i? 9c9s (r)] = E[^ (r)] - (i? 9c3s (r)) 2 = (24) 

Var[/\ s9c (r)] = E[i% 9c (r)] - (R 9s9c (r)) 2 = (25) 

Varfo^Cr)] = E[£ 2 s9 »] - (i? 5s9s (r)) 2 = Var^Cr)] (26) 

Var^Cr)] = E[i? 2 c9 >)] - (i? 9c5c (r)) 2 

_ (l + Jo(47r/ 1 T))(l + Jo(47r/ 2 T))+2Arj2(27r/ 1 r)(Jo(47r/ 2 r) + l)+2MJ 2 (27r/ 2 r)(J 

8NM 

+ "^iV 2 " 1(5 + 2M 2 ^ 2iV 2 M 2 l y A + ^ 1 ' T )w 2 ' T )J 

l + Jo(47r/ 2 r) + 2MJ 2 (27r/ 2 r) c ^ 1 + J (47r/ 1 r)+2jVJ 2 (27r/ 1 r) 
4i^M a/l ' T) iiVM 2 " ah ' T) (27) 

where = £n=i (E[A n (r)C n (r)} - E[I re (r)]E[C n (r)]) and Vg 6 = E!=i (E[5 m (r)/3 m (r)] - 

E[B m (r)]E[D m (r)}). 
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B. Second-order statistics for Simulator D 

Autocorrelation, cross-correlation, autocorrelation of the complex envelopes, and autocorrelation of the 
squared envelope for Simulator D are given below. Steps for the proof are presented in Appendix |B1 

Jo(27r/ir)J (27r/2r) + K cos (2nf 3 T cos(</> 3 ) 



Rh c h c { T 

R h c hA T 

R h s h s i T 
Rh s hA T 

Rhh{r 



E[h c (t + T)h c (t) 
E[h c (t + T)h s (t) 



2(1 + K) 
isT sin (27r/ 3 Tcos0 3 )) 



2(1 + K) 

■ E[h s (t + T )h s (t)] =R hchc {r) 

■ E[h s {t + T)h c (t)] = -R hchs (r) 

[h(t + r)h* (t)] = i (R hchc (r) + R hshs (r) - jR hchs (r) + jR hshc (r)) 

_ J (27r/ir) J (27r/ 2 r) + KeP^^^ 
1 + K 

R lhm 4r) = E[h 2 c (t)h 2 M + h 2 s (t)h 2 M + h 2 c (t)h 2 s ( Q ) + h 2 (t)h 2 ( 6 )} 

R\g\>\ g \>(T) + 8JTi? gege (r) cos (2tt/ 3 t cqs(0 3 )) +8i^+4K 2 
4(1 + If) 2 

For sufficient N and M, the expression of (l34l can be simplified as 



%| 2 H 2 ( r ) = JTTK) 2 ^ + J ° 2(27r/lT) + J o( 2 ^2T) + J 2 (2 7 r/ 1 r)J 2 (27r/ 2 r) 



The time-average correlations for Simulator D can be expressed as 

o , \ RgcgA T ) + -f^cos(27r/ 3 rcos(^ 3 
Rh c hA T ) = 



Rh c hA T ) = Rh s hA T ) 



Rh s h s (r) 



2(1 + JO 

X sin (27r/ 3 r cos(<; 
2(1 + A') 

A?»fl 5 ( T ) + K COS ( 27r /3T COS(0 3 )) 



2(1 + A") 

The variance of the time-average correlations for Simulator D can be derived as 



Var[/Wr)] = E{R 2 kk {t)] - (iWO) 



2 _ Var[fl gege (r)] 



4(1 + Kf 



V a r[R hshs (T)]=nR 2 hshs (T)} - (iWr)) 2 " Var fo-*^ 



4(1 + K) 2 



\ a t[R hohs (r)] = E[Rl hB (T)} ~ (Rh^)) 2 = 
Var[iWr)] = E[£ 2 (r)] - (i2h. h .(r)) 2 = 



(28) 

(29) 

(30) 
(31) 

(32) 

(33) 
(34) 



+ 2K(l + Jo(2vr/ 1 r)Jo(2^ 2 T) cos(2^/ 3 r cos(<£ 3 ))) + K 2 ). (35) 



(36) 
(37) 
(38) 

(39) 

(40) 

(41) 
(42) 
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IV. Simulations Results 



In this section, we present extensive Monte Carlo simulation results on the statistical properties of the 
proposed simulators. The simulation results are obtained using N s = 10 6 samples in a 10s time duration 
with sampling period T s = 10~ 5 s. The frequencies in all simulators are /i = /2 = WOHz, P = Q = 16 
in Simulators A and C, M = N = 16 in Simulators B and D, and <pi = 4>\2 = § in Simulators C and 
D. All plots are obtained using one trial unless stated otherwise. We also simulate the double -ring |[T3l 
and fixed-to-mobile channel models lfT4ll . ifTTl , |[27l to demonstrate the difference in statistical behavior 
of cascaded Rayleigh fading and single Rayleigh fading distributions, 

A. Probability density functions for channel envelopes 

The distribution of channel envelopes for cascaded Rayleigh channels without and with LOS has 
been discussed in lfl2l . ll28l . ||29l . Recognizing that the average power for simulators A and B is one, the 
probability density function (PDF), the cumulative density function (CDF) of the envelopes of Simulators 
A and B are given, respectively, by (major steps are provided in the Appendix O for convenience of 
readers): pz{z) = 2zKo(z\/2j, and Pz(Z < z) = 1 — \^2zKi(z^^, where z represents the envelope 
of Simulators A and B, i.e., \g(t)\, Kq and K\ are, respectively, the zero-order and first-order second 
kind modified Bessel function. Using (12) in | |29l with <t,- = , 1 . for i = 1,2, the PDF for the 
envelope of Simulators C and D is given by 



where z = \h\, Iq is the zero-order first kind modified Bessel function. 

Fig. ?? shows the simulation of the PDF of envelopes for Simulators A and B. The plots indicate 
that distributions of envelopes of both simulators match the theoretical curves very well. For comparison, 
simulations of the PDF for the double -ring M-M mobile channel simulator fT3| is plotted. While the 
envelope of the double -ring simulator follows Rayleigh distribution, Simulators A and B follow cascaded 
Rayleigh (worse than Rayleigh) distribution. 

Fig. ?? shows the simulation of the PDF of envelopes for Simulators C and D with different values of 
K. For the special case K = 0, Simulators C and D becomes, respectively, Simulators A and B up to a 
factor of This figure shows that the distribution of envelopes of both simulators match the theoretical 
curves in d43l ) for all values of K. 




VT+K 



VT+K 



(43) 
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B. Second-order statistics 

Autocorrelation and cross-correlation of Simulator B. The autocorrelation of Simulator B is shown 
in Fig. ?? using one and three trials. This figure suggests that Simulator B converges to the theoretical 
autocorrelation in one trial, and the difference between the simulation and the theoretical (fl3T ) is smaller 
in three trials. The cross-correlation is shown in Fig. ?? using one, three, and five trials. While simulated 
cross-correlation using one trial converges to the theoretical ( fTTT ) expression reasonably well, the simulated 
cross-correlation converges to the theoretical expression more closely when more trials are used. 

Autocorrelation and cross-correlation of Simulator D. The autocorrelation of Simulator D with 
different values of K is shown in Fig. ??. These plots indicate that the simulated autocorrelation converges 
to the theoretical expression very well in one trial. It can also be observed that the LOS component 
dominates the autocorrelation as r increases. This is because the Bessel function approaches zero when r 
is large, and the autocorrelation thus approaches K cos (2^/3 cos(^3)) (2(1+K)) . The cross-correlation 
for various values of K is shown in Fig. ??. Again, the plots indicate that the simulation results converge 
to theoretical expressions in (128T ) and (DTI ) in one trial. 

Fig. ?? plots the autocorrelation of the squared envelope for different values of K. The simulated 
results converges to its theoretical in (l34l ) in one trial. The plots suggests that the autocorrelation of the 
squared envelope approaches 1 when r is large, regardless of K. As r — > 00, all Bessel terms in (|34l ) 
become zero, which results in R\ h \2\ h \2(r) = K ^^ 2 ^t 1 = 1- 

Variance of the time-average correlations for Simulator D. The variance of time-average autocor- 
relation and cross-correlation for Simulator D with K = and 10 are obtained by averaging 200 samples 
for each r. As shown in Fig. ??, the variance time-average autocorrelation and cross-correlation are, 
respectively, in the magnitudes of 10 -4 or 10~ 5 . Since the plot for K = reflects a variance of the 
time-average correlations for Simulator B up to a factor of 2, the small variance in the figure confirms 
that all proposed simulators have satisfactory convergence performance. 

C. Higher-order statistics: LCR and AFD 

The level crossing rate and the average fade duration are two important statistical properties related to 
channel dynamics. At a specified level R, the LCR for a channel envelope is the rate (in crossings per 
second) at which the channel envelope crosses in the positive (or negative) direction ifTOl , ||27l , 11301 . The 
AFD is the average time duration that the envelope remains below the level R E71 . |[30l . Approximations 
for LCR and AFD are plotted in 1311 , ||32l for the cascaded Rayleigh fading channels. Here, we present 
plots by numerical evaluation of the exact LCR and AFD functions. 
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LCR for Simulator B. The LCR for the cascaded Rayleigh channel without LOS has the following 
form |29l, El: 

L\ g \(R) = (2nV2) *Rfi J™ x~ 2 Va 2 R 2 + x 4 e' 3 ^ dx (44) 

where a = j^, and R is the level. In the simulations of LCRs, the normalized LCRs and normalized signal 
level are used O, ifTOl . lfT3l . |[27l . For simulator B, the normalized LCRs is Ligi/j -1 , and normalized 
signal level is defined as p = -j=, where \fV is the root mean square (rms) envelope level for the 
channel, and V = 2. As shown in Fig. ??, the simulated LCR matches well with the theoretical (l44l . 
For comparison, we also simulated the LCR for a fixed-to-mobile channel lfl4ll . It can be observed that, 
for given signal levels, the LCR for Simulator B is more likely higher than that for fixed-to-mobile and 
double-ring models. This can be explained by the higher dynamics of the statistical properties in cascaded 
Rayleigh fading than a single Rayleigh fading distribution. 

LCR for Simulator D. The LCR for Simulator D is expressed as (29l, (331: 



3 /"CO r7T 

L W (R) = (l + K) 1 RfJ / x- 2 ^a? yi {R,d) 

JO J — 7T 



y 2 (R,#,x) \\ f2_ yi+g^+aia*!) 



x \ y 2 (R, x) [I + erf ^ ) ) + \j ^e~^^~ j ^ v 1+K {x+ ^^ ) dMx (45) 
where erf(-) is the error function 04l . and yi(R, ?9) and y2(R,$,x) are defined, respectively, as 



y 2 {R,#,x) = — - . (47) 

h^/VT+K^a? yi {R^) + x A 

The LCRs of Simulator D with different values of K are shown in Fig. ??. The normalized signal level 

is p = R, since the rms envelope level for Simulator D is one. While the plots indicate a good match 

between the simulated LCRs and the theoretical expression in d45l ), it is also observed that the higher 

the K, the lower the observed LCRs. 

AFD for Simulator B. Using (l44l and the CDF for \g\, the AFD for Simulator B can be obtained 

as El: 

T ]g] (R) = Pz ^ Z =^ = (l - V2RK X (Ry/2) ) f /iWWI x~ 2 Va 2 R 2 + x A e~ S ^dx^j (48) 

In the simulations, normalized AFD O, OH, US, O, H3 is plotted. For Simulator B, the normalized 
AFD is Tj/ji/i. The plots in Fig. ?? shows agreement between the simulated AFD with the theoretical (|48T ). 
For lower level of signals (p < dB), the AFD of Simulator B is longer than the AFDs of fixed-to-mobile 
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and double -ring models; and for higher level of signals (p > dB), the AFD of Simulator B is shorter 
than those of fixed-to-mobile and double-ring models. This indicates that cascaded Rayleigh fading is 
more severe fading than a single Rayleigh fading. 

AFD for Simulator D. The AFD can be calculated as T\ h \(R) = ^fjq , where L\ h \(R) is given 
in ( |45T ), and the CDF of the envelope z = \h\ is given by 



where Jo an d I\ are the modified Bessel function first kind, zero-order and first-order, respectively. Fig. ?? 
shows the AFD for Simulator D. For lower level of signals (p < dB), the value of AFD decreases 
when the values of K increase (stronger component of LOS), and for higher level of signals (p > dB), 
the value of AFD increases when the values of K decrease (weaker component of LOS). 

D. Complexity analysis 

In this section, we examine the complexity of our proposed channel models. Since simulators are 
sparse for M-M cascaded Rayleigh fading channels in the literature, we compare Simulators A and B 
with the double-ring simulation model, one of the well-known simulators for M-M Rayleigh fading 
channels. It is worth to mention that double-ring and the proposed simulators are for different categories 
of M-M channels, the comparison is solely for the illustration of complexity and performance of the 
proposed simulators. The number of additions required for Simulators A and B to generate one channel 
sample are, respectively, 2{P + Q) and 2(M + N), and 2MN additions are required for the double -ring 
model (same notations (M and N) are used in lfl"3ll ). Besides additions, one multiplication is required for 
both Simulators A and B, and no multiplication is required for the double-ring model. Table U lists the 
averaged CPU time elapsed when generating one channel sample for Simulators A, B, and double-ring. 
A HP Compaq 8510p computer, Intel (R) Core (TM) 2 Duo CPU T7500 @ 2.20GHz, is used with 
M = N = P = Q = 8 for all simulators. It can be seen that the elapsed time to generate one sample 
for Simulators A and B is much shorter than that for the double -ring model. A relative CPU time is also 
listed, with To being the referenced time for the double-ring model. 

To investigate the performance of convergence, we examine the mean square error (MSE) between 
simulated autocorrelation and the theoretical for Simulators A, B, and double -ring model. These simulators 
have an identical theoretical autocorrelation expression (fT3T ). The MSEs are obtained using one trial for 
different complexity levels. As shown in the first four rows in Table JIJ if the same number of additions 
are used to generate one channel sample, then both Simulators A and B have much better convergence 
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performance with significantly smaller MSE than the double-ring model, and Simulator B provides the 
best convergence performance. If the same number of sinusoids are used (in this case, Simulators A and 
B have much less complexity), the convergence performance of Simulators A and B is still comparable 
to that of the double-ring model. Fig. ?? plots the autocorrelations of Simulators A, B, and double -ring 
model for the scenario listed in the second row of Table [TTJ whereby the number of additions to generate 
one channel sample is 200 for all simulators. The plots indicate that Simulator B has a faster convergence 
rate than the double -ring model. 

TABLE I 

CPU TIME ELAPSED FOR ONE CHANNEL SAMPLE 



Model 


Computations for one channel sample 


CPU time 


Relative 
CPU time 


Simulator A 


2(Q + P) additions, 1 multiplication 


3.8s 


0.2 T D 


Simulator B 


2(M + TV) additions, 1 multiplication 


4.3s 


0.23 T D 


Double-Ring 


2MN additions 


18.8s 


T d 



TABLE II 

MSE AND COMPLEXITY 



Additions for One 
Channel Sample 


Simulator A 


Simulator B 


Double-Ring 


128 


5.34 x 10" 4 
{Q = P = 32} 


2.12 x 10" 4 
{M = N — 32} 


8.1 x 10" 3 
{M = N = 8} 


200 


3.19 x 10" 4 
{Q = P = 50} 


1.80 x 10" 4 
{M = TV = 50} 


4.69 x 10" 3 
{M — N — 10} 


288 


1.92 x 10~ 4 
{Q = P = 72} 


1.64 x 10 -4 
{M — N — 72} 


2.91 x 10" 3 
{M — N — 12} 


392 


1.78 x 10~ 4 
{Q = P = 98} 


1.28 x 10~ 4 
{M — N — 98} 


2.40 x 10~ 3 
{M — N — 14} 


32: Simulators A and B 
128: Double-Ring 


7.3 x 10" 3 
{Q = P = 8} 


2.8 x 10" 3 
{M = TV = 8} 


8.1 x 10" 3 
{M = TV = 8} 



V. Conclusion 

We have proposed statistical simulators for mobile-to-mobile channels, whereby the received signals 
experience cascaded Rayleigh fading with or without LOS. The simulators contain two individual 
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summations and are therefore easy to implement with lower complexity. Furthermore, the simulators 
provide faster convergence to all the desired statistical properties, including the pdf, autocorrelations, 
LCRs, and AFDs. Theoretical derivation of the time-averaged statistical properties and the corresponding 
variance are derived to confirm that the proposed simulators have good convergence performance. 
Extensive Monte Carlo simulation results on various statistical properties and complexity analysis are 
provided to validate the proposed simulators. While measurements and tests in various highly dense 
scattering environments confirm that the M-M channels may undergo cascaded Rayleigh fading (more 
severe than single Rayleigh fading), our proposed simulators can be used to simulate the underlying 
channels and reveal the corresponding statical properties. 
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Appendix 

A. Proof of (1151 ) and (1171 ). autocorrelation and cross-correlation for Simulator B 

From (O, the real and imaginary parts of Simulator B are given, respectively, by 

r-^— N,M 

g c (t)=Re(g(t)) = J— E (A n (t)B m (t)-C n (t)D m (t)) (50) 

n,m=l 



N,M 



g,{t) = hn(g(t)) = \/^rr E (A»(t)A»(t) + C n {t)B m (t)) (51) 



NM 

n,m=l 



For brevity of notations, we replace t + r with g in the proceeding appendices. The autocorrelation can 
be obtained by 



„ r N,M N,M 

E[g c (t)g c ( Q )] = —E £ (A n (t)B m (t) - C n (t)D m (t)) ]T (A^B^g) - C^D^g 

.n,m=l k,j=l 



N,N M,M N,N M,M 

E[ MQMe) E ^m(t) J B i ( e )+ E C n (t)C k (g) E D m (t)Dj(g) 



NM 

n,k=l m,j=l n,k=l m,j=l 

N,N M,M N,N M,M 

- E Mt)Ck(g) E B m {t)D 3 {g) - E C n (t)A k (g) E D m {t)Bj{ e )] (52) 

n,k=l m,j=l n,k=l m,j=l 

The cross-correlation is evaluated as 

„ N,M N,M 

E[g c (t)g s (g)] =E[— E (A n (t)B m (t) - C n (t)D m (t)) E (A j (g)D k (g) + C^B^g))] (53) 

n,m=l j,k=l 
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Since n , <D m , VE^ E [— ir,ir) are statistically independent and uniformly distributed for 
all n and m, we have E[E^ lin#fc A n (t)A k (g)] = EE^i,^ C„(t)C fc ( e )] = 0, and 
EES^^^B^e)] = E[E^ Wj ^m(*)^i(e)] = 0. It is easy to justify E [<7 c (*)<?c(f?)] = 
Jo(27r/ 1 r)Jo(27r/ 2 r) and E[g c (t)g s (g)] =0. 

Using |[34l (p. 420-421), the following identities are listed for the convenience of the proceeding proof 



of the simulators' statistical properties. 

N,N N 



1 N N 
E[ 2 4(i)^(e)]=EEA»(i)A,( ? )] =;K[E4(r)] =-^(2^) 



n,k=l 
N.N 



it 

N 



n=l 
N 



E[ ^ CnW^(^)] =EEC„(t)C n (e)] =^EEC'n(r)] ^Jo^r) 

n,k=l n n=l 

M,M M , Af ^ 

£ HB rn (t)5 J -(e)] = X; E[4p m (?)] = -EE B to (t)] =—J (2irf 2 T) 

m=X,j m=l m=l 

£ HD m (t)D i (e)] = X; nD m (t)D m (g)] = -J2 n£> m (r)] =— J (27r/ 2 r) 



m=lj 
N,N 



m=l 

MM 



m=l 



J2 E[An(t)C k (g)] = nB m (t)D 3 (g)\ = 0. 

mj=l 



n,fc=l 



B. Proof of (1281) ancf (1311) . autocorrelation and cross-correlation for Simulator D 
The real and imaginary parts of Simulator D, respectively, are 

h c (t) = Re(h(t)) = -j=t=Ulg c {t) + VKcos (L(t))\ 
h s (t) = lm(h(t)) = -^t=U^g s (t) + ^sin (L(t)) 



where L(t) = 2tt f^t cos(<ps) + 

1 



The autocorrelation is calculated as 



E[h c (t)h c (g) 



l + K 



;E 



\g c {t) + VKcos (L(t))) (\\gc{Q) + ^cos (l(q)) 



The cross-correlation of Simulator D can be obtained as 



E[h c (t)h s (t+r)]= T -^E ^I 5c (i) + ^co S (L(t))^y^(t+r) + ^ S in(L(t+r) 



(54) 
(55) 
(56) 
(57) 
(58) 



(59) 
(60) 

(61) 
(62) 



Notice that the phase <po in L(t) and L(g) is independent of other random variables in g c (t) and g s (t). 
Taking the expectation with respect to <fto and using the results in the autocorrelation and cross-correlation 
for Simulator B, one can obtain the autocorrelation and cross-correlation as specified in (|28T ) and d3TT >. 
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C. Proof of d 19b .' Squared envelope correlation for Simulator B 

The squared envelope correlation for Simulator B can be written as i?| g |2| g |2(r) = E[g 2 (t)g 2 (g)~j + 
E[g 2 (t)g 2 (g)] + E[g 2 {t)g 2 (g)] + E[g 2 (t)g 2 (g)] . The first term is expressed as 



E[g 2 (t)g 2 ( g ) 



N 2 M 2 



r N,M N,M 

E Y (A n (t)B m (t) - C n (t)D m {t)) Y (A u (t)B p (t) - C u (t)D p (t)) 

-n,m=l u,p=l 
N,M N,M 1 

x Y (Mq) B Aq) ~ C g (g)D r (g)) Y {A s (g)B J (e)-C s {g)D J (g)) . (63) 

q,r=l s,j=l 

Expanding (l63l and taking expectation with respect to the phases, we obtain E^(t)^(^)] = 
^(TaTb + T C T D + M^f. + 4k), where T x = E[y£ X n (t) Z% X v (t) X q {g) £f X s (g)], 
where X = {A, B, C, D}, and 

N,N N,N M,M M,M 

K = E[YMt)As{Q)Y C ^ t ) C ^e) Y B m(t)Bj(g) Y D P (t)D r (g)] 



(64) 



u.q 



Similarly, we have E[g 2 (t)g 2 (g)] = E[g 2 (t)g 2 (g) 
K[9 2 s(t)g 2 (8) 



4 

N 2 M 2 



p,r=l 



( Mi (T A + T C ) + ^ (T fl + T D ) - Ak) , and 



-(T j4 T D + T c T,3 + 



KPN 2 



N2M2\± a±l> -r o- j-±( n § h 4re). The term T A is evaluated as 

N N N,N N,N 

T A = E[YA 2 p (t)YA 2 J (g)]+E[ y Mt)Mt) E Mq)Mq) 

n,u,n^u q,s,q^s 



(65) 



P 3 

The first term in d65l ) is obtained as 

AT 



P 3 

The second term in (l65l) contains the following seven cases: 

• Case 1. n ^ u, q ^ s,n = q,u ^ s; 

• Case 2. n ^ u,q ^ s,n = s,q ^ u; 

• Case 3. n 7^ u,q 7^ s,u = q,n 7^ s; 

• Case 4. n ^ u, q ^ s,u = s,n ^ q; 

• Case 5. n ^ u,q ^ s,n ^ q,u ^ s; 

• Case 6. n ^= u,q ^ s,n = q,u = s; 

• Case 7. n 7= u, q 7= s, n = s, it = q. 

The value of e[^^ A n {t)A n {g) J^^s^n A s {t)A s (g)] is zero for Cases 1 to 5 and identical for Cases 



6 and 7, and E[^ A n (t)A n (g)j:^ n A s (t)A s (g) 



N 2 J 2 (2Kf lT )-J2: (E[A„(r 



case 6 or 7 



-. Therefore, 



'+NJ (4nf 1 T)+4N 2 J 2 (2nf 1 r) _ Z" ( E [^"( T )l) 
8 

obtained following the steps in (48) Appendix I |[35l . 



we have T A = ^-t^.w. T m»- Jo - WlT , _ ^„ v -,-.-, v _ ^ fect> ^ valug rf T ^ can ^ bg 
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o- l l u 'V 2N 2 +NJ Unf 1 T)+4N 2 j2(2nf 1 T) E» ( E Ml ) T . ■ . • Utf , . . . f 

Similarly, we nave Tc = — — iL -g J — — - J — L 2 — • It is straightlorward to justify 

E[c , „(t)] = E[CN-n{r)],n = 1, • • • , JV by the following steps: 

E[C„(t)] = f cos (27T/ir cos( 2n7r 7 J + ^ )) i-# 



4iV 



2vr 



2iV 



2iV 



COS 



/ \ 4iV 

(2vr/ 1 rcos(0))— tty 



AN 



cos (27r/irsin(^))— = e[CV_„(t)] 



(67) 



Denote f(/ 1>r ) = £n (E[i n (r)]) = En (Ep„(r)]) . Then, we have 

2N 2 + NJ (4nfiT) + 4iV 2 J 2 (27r/ir) e(/i,r) 



T A = T C 



T B = T D 



2M 2 + MJo(4vr/ 2 r) + 4M 2 J 2 (2vr/ 2 r) £(/ 2 , r) 



(68) 
(69) 



££(E[ 



Em ( e [Ai(t)]) - Inserting the results for 



where £(/ 2 ,r) 

E[ 5 2 (i)<? 2 (£»)], E[g 2 c (t)g 2 s (g)],E[g 2 s (t)g 2 c (g)], and E[(£(t)^( e )] into % P | 9 p(r), we obtain 
%l 2 l9l 2 ( r ) 



4 / N 2 \{ M 2 

WW, T A + To + - T, + T D + - 



A^+iV 2 J 2 (27r/ 1 r) + 



iVJ (47r/ir) 



Expanding (fTO]) yields ( fl9l ). 



^(/ 1 ,r))(M 2 +M 2 J 2 (27r/ 2 r)- 



MJ (47r/ 2 r) 



-£(/2,T 

(70) 



D. Proof of (134b : Squared envelope correlation for Simulator D 

The squared envelope correlation for Simulator D contains four terms as indicated in (133T ). The first 



term is evaluated as 



E[7£( t )fc2 (e) ] = 1 



if) 2 



(1 + K) 2 



E 



i(/ c (t) + v^cos J (J^9c(q) + \FK cos (L(e) 

^ 2 (i) + V2Kg c (t) cos (L(t)) + Kcos 2 (i(t)) 



-g 2 c {g) + V2Kg c {g) cos + K cos 2 



(71) 



Note that E[g c (t)} = E[g c (g)] = 0, E[cos (.£(#))] = E[cos (iL (g))] = 0,£ = 1,2, and E[ 5 2 (£)] 
E^fe)] =i2 flcffc (0) = l. We have 



Ep£(ty£(<?) 



^E[ 5 2(t) 5 2( e )] + 2K^ c (r)E[cos (L(t)) cos (L( e )) 



(1 + iT) 2 U 
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+ ^E[cos 2 (L(0) +cos 2 (L(£))] + A 2 E[cos 2 (L(t)) cos 2 (L(q))]J (72) 

The autocorrelation of the squared quadrature component E[/i 2 (t)/i 2 (^)] can be evaluated fol- 
lowing similar steps in obtaining E[/i 2 (t)/i 2 (£>)]• Using identities E^cos (-k(i)J cos (L(g))] = 
icos(2vr/ 3 rcos(^3)),E[cos 2 (L(i))] = E[cos 2 (L(p))] = ±, and E[cos 2 (l(t)) cos 2 = 
\ + | cos (47r/3T cos(03)), we summarize the auto- and cross-correlations, respectively, as 

xnr,2, „ 2 / n1 2E[ ff 2 ^) ff 2 (£)]+8A^R gx9 >)cos(^^^^ 

E[^(t)^(^)j = 8(1 + *T) 2 

(73) 

wk2 , u2/ Al w r,2^z,2/ M 2E[ 5 2( t ) 5 2( e )] +4A+2A 2 -A 2 cos(4vr/ 3 rcos(03)) 



8(1 + A) 2 

where x = {c, s}. Inserting d73l and (l74l into (l33l yields (l34l . 



£. Proof of (127b : Variance of time-average correlations for Simulator B 

The variance of the time-average autocorrelation of the real part for Simulator B is Var[i? 9c9c (r)] 
E[Rg cgc (r)] - (Rg c gXr)) 2 - The first term is evaluated as 



1 



N,M 



NM 



nR'l 9 M= J^2 E [ Y,{ A n^)Bm(r) + C ri (r)D m (r)) Y,{ A p(r)B q (T) + C p (r)D q (r)) 



1 



n,m=l 
N,N 



4N 2 M 2 



E 



MM 



p,q=l 



+2E 



N,N 



MM 



+ E 



Ln,<j=l m,p=l 
N,N MM 

£ C n {T)C q (r) J2 D m (r)D p (r 

-ra, q=l m,p=l 

We have the following identities: 

N,N N 

E[ Mr)A q {r)\ =iY 2 J 2 (2^/ 1 T)+E[^i n (r)i n (T)] -^E[i n (r)]E[I n (r 



Ln,q=l m,p=l 

(75) 



N 



n, q =l 



71=1 



71 = 1 



Joi^Jir) + — H ?(/i,r) 



E[ £ C B (r)C,(r)] =E[ X i„(r)A,(r 



7i,g=l 
MM 



n,q=l 
MM 



E[ X B m (r)^(r)]=E[ £ 4 (r)D p (r )] = M 2 J 2 (2vr /2 r ) + M + M J ° ^ 



7Tl,p=l 

N,N 



m,p=l 



E[ X 4(r)C 9 (r)] =iV 2 J 2 (27r/ 1 r) + y 

71,5=1 



AC 



(76) 
(77) 

e(/ 2 ,r) (78) 
(79) 
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M,M 



E[ B m (T)D p (T)]=M 2 J$(2irf 2 T) + V Sf) 

m,p=l 



(80) 



where V Ad = £n=i (E[^ n (r)C n (r)] - E[A»(r)]E[C n (r)]), = Em=i (E[J3 ro (r)D m (r)] 

E[S m (r)]E[£> m (r)]). Inserting |76]> to dHHl) into |75]), we obtain 

x (M 2 J 2 (2tt/ 2 t) + f + ^iMl^ _ e(/ 2 ,r))+ 2 (iV2j 2 (2vr/ 1 r) + ^)(M 2 J 2 (2vr/ 2 r)+^ 

= ^^((27V 2 J 2 (2^/ 1 T)+iV+iVJo(47r/ 1 r)-2 £(/ 1} t)) 

x(2M 2 J 2 (27r/ 2 T)+M+MJ (47r/ 2 T)^ 

(81) 

Recall (R gc9a (r)) 2 = J$ (2n fcr) J$ (2n f 2 r) . We have 
1 



Reorganizing (1821 . one can obtain 

The variance of time-average autocorrelation of the imaginary part of Simulator B is Var[12 gsffo (r)] 
E[^ 2 Ss (r)] - (i? 9s9s (r)) 2 . The first term is evaluated as 

N,M N.M 



Var[i? 9c9c (r)] = ^^( 2 M 2 J 2 (2vr/ 2 r)(iV + iVJ (47r/ir) - 2 ^r) + 2V^) 
+ 2iV 2 J 2 (27t/ 1 t)(M + MJ (47T/ 2 r) - 2 £(/ 2 , r) + 2Vg 6 ) 

+ (N+NJo^irfxT) - 2 e(/i,r))(M+MJ (4vr/ 2 r) - 2 £(/ 2 , r)) +4V^Vg 5 ) (82) 



4M 2 N 2 
1 

4iV 2 M 2 



n,m=l Pi9=l 

AT, TV MM r AT, AT M,M 



E 



Y,A n {T)A q {r)Y, i) rn{T)D P {r) 

L n,g=l m,p=l 
AT, AT M.M 



+2E 



^I n (r)<5 ? (r)^5 m (r)D p (r) 

-n,q=l m,p=l 



+ E 



Y, C n (r)C q (r) Y B m (r)B p {r) 

-n,q=l m,p=l 



(83) 



We have Var[4, 9s (r)] = Var[i? 3c9c (r)]. 



Proof of d39b ." Variance of time-average correlations for Simulator D 

We present the derivation of the variance of the time-average autocorrelation of the in-phase component 
in d39l ), while the same steps can be applied to obtain the variance of the time-average autocorrelation of 
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the quadrature component. The variance of the time-average autocorrelation of the in-phase component 
for simulator D is Var[i?^ c /j c (r)] = E[i?| c/lc (r)] — (Rh c h c { T ))' The first term is evaluated as 

E fe,( T )i^[( ^ (T)+ X ( ?) /3TC ° S " 3)) ) 2 ' 

= 4(1 ^ )2 (E[fig cgc (r)] + 2A-COS (27r/ 3 rcos(03))E[fl ( , el , c (r)] + (K cos(2^/ 3 r cos(</> 3 ))) 2 ) (84) 

Recall E[£ ffofle (r)] = J (27r/ir) J (27r/ 2 r). Using E [fl^ ( r )] obtained in Appendix E and R hc hS T ) 
from Appendix |B] it is easy to obtain 



G. Proof for PDF of envelope for Simulators A and B 

The pdf and phase distribution of Simulators A and B are readily available in the literature |[28l . Denote 
9 = 9192 = 9c + j9s, 91 = xi +jx 2 , 92 = Vi + jV2, %i,x 2 ~ Af(0, of), and y 1 ,y 2 ~ Af(0,of) are i.i.d. 
zero mean Gaussian normally distributed. The joint PDF for the real and imaginary parts of g is given 
as Pg c g a (g c , 9s) = 2irl'(a A ^° ( ^ai^ ''^) ' ^° ^ nc ^ ^ e P^ °^ ^ enve l°P e > we transform the Cartesian 
coordinates (g c ,gs) to polar coordinates (z,7]), where z = y^|~+~gf, and r/ = arctan^l^^). The 
resulted transformation Jacobian is 1/z. The joint pdf of the envelope and phase is given by p z ,-q{z, if) = 
z Pg c gA9c,9s) = 2™yt K °C^k\ The PDF of the envelo P e is obtained as p z (z) = -^K {^f^j. 
Recognizing c?i = a 2 = l/yv^ for Simulators A and B, one can obtain the PDF of their envelopes. 
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